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Abstract 

We derive the off-shell photon propagator and fermion-photon vertex at one-loop level 
in Maxwell-Chern-Simons quantum electrodynamics in arbitrary covariant gauge, using four- 
component spinors with parity-even and parity-odd mass terms for both fermions and photons. 
We present our results using a basis of two, three and four point integrals, some of them not known 
previously in the literature. These integrals are evaluated in arbitrary space-time dimensions so 
that we reproduce results derived earlier under certain limits. 

PACS numbers: 11.10.Kk,11.15.Bt,11.16.Yc 
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I. INTRODUCTION 



When local gauge symmetry is what drives theories of fermion and boson interactions, 
mass terms for the latter are forbidden. Thus, the mass of these particles should be explained 
with different arguments such as the spontaneous breaking of symmetry through a Higgs 
mechanism. The recent discovery of a boson with 125 GeV mass and quantum numbers 
similar to those of the Higgs boson [l| brings this mechanism to a physical reality in the 
high-energy realm. Nonetheless, if leaving local gauge symmetry untouched is desired and 
yet one wishes to allow the gauge boson to acquire a mass, one could rely on topological 
quantum field theories such as Chern-Simons (CS) theory [2]. The CS theory has been 
put forward in (2+l)-space-time dimensions in several contexts of condensed matter physics 
(see, for instance, Ref. jsl and references therein), quantum gravity {4] and string theory jsl. 
The Lagrangian for such a term, although not manifestly gauge invariant, possesses a gauge 
invariant action. Particularly, considered along with the Maxwell theory for photons, the CS 
term provides a topological mass for these gauge bosons, which preserves gauge invariance. 
Moreover by adding matter fields, the theory becomes richer in the sense that some features 
are apparently unique to (2+l)-dimensional theories, such as anyon excitations, generalized 
parity, and so on. 

In this paper, we study the quantum electrodynamics of four-component fermions and 
photons in two space and one time dimensions coupled to a CS term, QED3. The action of 
this theory is gauge invariant, though parity and time reversal might be explicitly broken. 
Fermion and boson mass terms play an intricate role, because in addition to the usual fermion 
mass term, a parity violating mass term can be included. Such a term would radiatively 
generate a CS term and vice versa (see, for instance, Refs. Specific cases of the 

underlying Lagrangian have been used in the description of several condensed matter physics 
systems, once fermions have been integrated out, such as high-T^ superconductors , the 



quantum Ha! 



fsee Ref. 



1 effect jsl, [lo| and, more recently, graphene Uj and topological insulators 



12| for a recent review). In high-energy physics, the use of QED3 has been 
connected to the study of dynamical chiral symmetry breaking and confinement, where it 
provides a popular battleground for lattice and continuum studies 13|. In particular, the 
CS term allows the possibility of chiral and parity symmetry breaking. Within the two- 
component fermion formalism, it was observed that while chiral symmetry can be broken. 
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parity symmetry remains untouched [ij]. Moreover, this term induces a detaching of chiral 



symmetry breaking-restoration and confinement-deconfinement transitions 



15|, leading to 



a phase resembhng quarkonia, where chiral symmetry is restored, but confinement takes 



place 



16| . Confinement and screening effects have been studied in this theory [l7| through 



the behavior of the vacuum polarization. 

A key issue in these t ype of studies is to address the gauge covariance properties of 



Green functions. In Ref. 



18|, the gauge covariance of the fermion propagator was estab- 



lished through the use of the Landau-Khalatnikov-Fradkin transformations 19|. In the same 



spirit, in Ref. 



20|, the fermion and photon propagator, and particulary the on-shell fermion- 



photon vertex were calculated and some limits explicitly considered. The same three-point 



function was studied in Refs. 



21 



22| . the later study carried out in the light cone gauge. 



Interestingly, in all these works, the corresponding anomalous magnetic moment for elec- 
trons is identified directly. Here we extend these findings computing the off-shell forms of 
these Green functions in arbitrary covariant gauges. The seminal work of Davydychev et 
al. 23| on the one-loop quark gluon vertex already provides a guide towards the master 
integrals involved in the calculation of the two- and three-point Green functions. Moreover, 



in Refs. 



24 



25| . the fermion-photon vertex for the parity preserving version of QED3 was 



calculated for massless and massive fermions, respectively. Reference 25|, in fact, provides 
the most general form of the non-perturbative fermion-photon vertex in QED3 guided by 
perturbation theory. In this work, we attempt to generalize these findings incorporating an 
explicit topological photon mass. In this task, we are faced with new master integrals which 
were not previously reported in literature, as far as we know. These involve an explicit 
gauge boson mass, and therefore we evaluate them in arbitrary space-time dimensions, so 
that these can be useful in other contexts. For us, these steps are fundamental in order to 
pursue further analysis of certain physical observables sensitive to the infiuence of the CS 
term, as will be reported elsewhere. 

The paper is organized as follows: In Sec. II, we consider four-component spinors to 
describe the QED3 Lagrangian with parity conserving and violating mass terms for fermions 
and photons. In Sec. Ill, we calculate the photon propagator at one-loop level for the 
most general Lagrangian. The one-loop correction to the fermion propagator is discussed in 
Sec. IV. In Sec. V, we derive the fermion-photon vertex at one-loop level under a general 
decomposition of the vertex we propose. Given this structure of the vertex, new master 
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integrals arise that we calculate and study some of their limiting cases, including the massless 
case for both fermions and photons. At the end, we summarize our findings in Sec. VI. 
Details on the calculation of the new master integrals, are presented in the Appendix. 



II. MAXWELL-CHERN-SIMONS QED3 

Let us start by presenting our model and conventions. We work with four- component 
spinors and thus with a 4 x 4 representation for the Dirac matrices. We choose to work 
in Euclidean space, where the Dirac matrices satisfy the Clifford algebra {'jpy'ju} = 25^j,. 
Once we have selected a set of matrices to write the Dirac equation, say {70,71,72}, two 
anticommuting gamma matrices, namely, 73 and 75, remain unused, leading us to define two 
independent global phase transformations of the chiral-type for the spinors: tl' — e'^'^'^^tp and 
ip — >■ e^^'^^ip , where (3 and a are arbitrary real numbers. Consequently, there exist two types 
of mass terms for fermions that can be considered in the Lagrangian, the ordinary mass 
term rrieipip and a new mass term nioipTip with r = ^[73,75], sometimes referred to as the 



Haldane mass term [26[ in condensed matter physics. The former violates chirality, whilst 
the latter is invariant under "chiral" transformations. Moreover, meipip is parity invariant 
but TrtoipTil) is not. This allows parity and time reversal to be broken in our theory, which 
is described by the Lagrangian 

- 11^ 
£ = ^(i ^ + e 4 - me - rm^)^ - -F^^F^^ - —{dpApf + -e^^pApF^p, (1) 

where the first term corresponds to the Dirac Lagrangian with parity-even and -odd mass 
terms and fermion-photon interaction; the next one is the Maxwell term followed by the 
covariant gauge fixing term, ^ denoting the gauge parameter. The last term is precisely 
the CS term, with 6 being the topological mass for the photons. Rather than working with 
parity eigenstates, we find it convenient to work in a chiral basis. For this purpose, we 
introduce the chiral projectors 

X± = lil±r), (2) 

which have the properties 

X± = X± , X+X~ = , x+ + X- = 1 . (3) 



4 




FIG. 1: One-loop correction to the photon propagator. 



Then, the right-handed tp+ and left-handed tp- fermion fields are ip- 



± 



fashion that the chiral decomposition of the fermion propagator becomes 18 



X±4', in such a 



SP(P) = -4^X. - , (4) 



with m± = rrif. ± rrio. Furthermore, we identify the photon propagator associated with the 
Lagrangian ([T]) as 

One can identify from the transverse piece of the propagator the role of ^ as a photon 
mass. It also plays a role in the parity-odd piece (proportional to e^^p), but it is absent in 
the longitudinal part. In the following section we calculate the one-loop correction to this 
propagator. 

III. THE ONE-LOOP PHOTON PROPAGATOR 

The one-loop correction to the photon propagator is depicted in Fig. [1] Using Feynman 
rules, it corresponds to the expression 

n^vlp) = J -^^SP{w)e-ff,SP{w - p)e-f^ , (6) 

which we define in terms in the left- and right-hand side projections as 

U^u{p) = n(+) ^^{p) + n(_) ^^(p) , (7) 

where 

(±) pAP) - J (2 7r)3 [w^ + ml] [{w - pf + ml] ' 
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To evaluate these expressions, we propose the following decomposition: 

a 



n 



(±) fJ.U 



ip) 



27r2 



where we have written the propagator in terms of 



K 





Pf - 


\- mW 










Pf- 










[w"^ + mj_][(M; — 


pY- 





(9) 



These can be evaluated straightforwardly, using standard techniques (see the Appendix for 



further details and also Ref. |27t]). Following Ref. [l8 

/(p, m 

so that the basis integrals may be expressed as 
K 



arctan 
p m 



we define 



P 





71 1 






2 








. 2 V 



+ vr 



2PjJ^ 
p2 



+ 4m^ 



(10) 

(11) 
(12) 

(13) 



Using the results (]TT|) -( 1T3|) into the decomposition ([8]) and making the replacement m± m 
we obtain the photon propagator in accordance to the results in Refs. [3, H] • This completes 
our calculation of the one-loop photon propagator. Next we calculate the fermion propagator 
at one-loop. 



IV. ONE-LOOP FERMION PROPAGATOR 



The one-loop correction to the fermion propagator has been previously considered by some 



of us in Ref. 18[ , and corresponds to the diagram in Fig. |2l To calculate this correction, we 
write the one-loop fermion propagator as 

S^P{p) = S^l\p) + S^l\p), 
6 



(14) 




FIG. 2: One-loop correction to the fermion propagator. 



where 



The functions A^^\p) and Bj^'{p) are thus obtained as 



S±{P) 



A'1\p) p^+ B'i\p) 



(15) 



2'Ka f <P 



w 



Tr 



Sl\w)^,m{w-p-Ox± 



'^,Si\w)^M^}{w-p-i)x± 



(16) 



Performing the traces, after suitable manipulations, these functions can be expressed in the 
form (compare against Eq (26) of Ref. 18| ) 

{p^ — rn?)'^ 



Al\p) 



[X2(/,m2)-X(0,m|,/)] 



a 



{p^ -m^ fl{0,ml,p'^) - [{p 



2'ir^p'^ 



{p' + m'- e')I{e\ mi,p') - {p' + m')I{0, m\p') - T{e') 



Bl\p) = ^X(0,ml,/) + ^X(e^m|,p2) 



27r2 
a 



271^6 



{p' + e' + mi)x{e\ mi,p') - {p' + m^)x(o, m',p') - r(r ) 



(17) 



where the master integrals are 



r(m|) 



X{e\ml,p') 
X2(m^,p^) 



9 I 2 ' 



(Tw 



(w'^ + 6^2] [(p — wY + mj_] ' 
1 

w'^[{p — wY + 



(18) 
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These integrals are performed using standard techniques (see Appendix), and give 

^^('"-^ > = -wtA^ ■ '''' 

Inserting the above results into the decompositions f|T71) and f|T5l) . we obtain the one- loop 



correction to the fermion propagator, in accordance with Eq. (27) of Ref. 18|. This com- 
pletes the one-loop correction of the two-point Green's functions corresponding to the La- 
grangian ([T]). Next we calculate the one-loop correction to the fermion-photon vertex. 

V. ONE-LOOP FERMION-PHOTON VERTEX 

At C(q:), the fermion-photon vertex is depicted in Fig. |3l and corresponds to 

r^(A;,p)=7M+A/.- (20) 
Using the Feynman rules, A^ is simply given by 



(21) 



Given the properties of the projectors in ([3]), the crossed structures vanish and therefore we 
have only integrals with the same mass propagators (m+ or mJ). Given this, we can write 



A^ as 



A;. = A(+) ^ + A(_) ^ (22) 




FIG. 3: One-loop correction to the fermion-photon vertex. 
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where 

A(±) t, = ^{ [la^lf, }ila + m±{4kf, + Ap^ - ^7^ - 7^/{f) - m^7j J^^\e, m±) 

H[lf.K^'\rn±) - [7,;^7^ + 7^^7, - 2m±5^,] Ji^)("^±) 
+ [7;.^7m^7a + "^±(7i^^7m7a + 7;^7/. /i^7A) + "^|7;.7m7a] /i?("^±)] 
-7^[irW(^±) - ^V(°)(^,^±)] - [7.^7m + 7m^7. - 2m±5^,] Ji^^^.^i) 
+ [7i^^7M/'^7A + "^±(7i'^7m7a + 7i^7/^ /'^7a) + m|7 7,7a] 

-^ez.Ap[- [7z^^7m 7p 7a + 7i^ 7p 7m /^7a + ^±(7^ 7^ 7^ 7a + 7^. 7m 7p 7a)] J^^^ {0, m±) 
+7;.7<7 7m7p7a4^H^>"^±) + [7;^^7m/^7a + "2±(7^^7^7a + 7i^7p /^7a) 
+m|7.7p7A]4'n^,"^±)]}x±- (23) 

In this expression we have introduced the basis integrals K^'^\ J^i^"*, J^?/\ 



IfiJ , Lfj, , L)j,l which we define as 



(2) r(2) 

1 



[{p — wy + [{k — wY + rrtj-] 
1 



j(i)(m±) = /" c/3y; 



J^^)(e,m±) = j d^w 
J!.'jim±) = [ d'w 



^'j(^,m±) = j d^w 
/(0)(m±) = [ d^w 



ll^\m±) = j d^w 
= [ d'w 



w"^ [{p — wY + 


[(fc — wY 
1 


+ ml.] 




[w^ + e^] [{p - 


U7)2 + m^] [(A; 


— wY - 




w"^ [{p — wY + 


m^] [(A; — wY 


+ ml.] 




[^2 + ^2] [(p _ 


wY + fn'^] [{k 


— wY - 




[(p — wY + 


m^] [{k — wY 

W^Wu 


+ ml.] 




[^2 + ^2] _ 


wY + m^] [{k 
1 


— wY - 




[(p — wY + 


[{k — wY 


+ 




[(p — wY + 


m^] [(A; — wY 


+ ml.] 




[(p — wY + 


m^] [{k — wY 


+ 


5 



L^^V^ m^) = f d^w ^ 

^ ' J W^W^ + ^2] [(p - w)2 + [(A; - u;)2 + ' 
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41H^,^±) = / d'w W—^- (24) 

^ J [w^ + [{p — wY + mj_J [[k — wY + m^J 

From this list of integrals, those which depend only on fermion masses have been cal- 
culated in {25! and for them we just replace m — m± to obtain the results needed in the 
calculation at hand. The rest correspond to two-, three- and four-point integrals that de- 
pend both on a fermion and a gauge boson mass, therefore they might be useful in theories 
in other space-time dimensions. We solve these in arbitrary dimensions in the Appendix. 
Also, it is important to mention that these results are crucial to our interest in the impact 
of the CS term on certain physical observables related to the Lagrangian ([1]). 

These results along with the decomposition proposed in Eq. complete the calculation 
of the one-loop correction of the fermion-photon vertex. We summarize our findings in the 
next section. 

VI. SUMMARY AND CONCLUSIONS 

We have derived the one-loop correction to the photon propagator, fermion propagator 
and the fermion-photon vertex in Maxwell-Chern-Simons in QED3. In doing so, we have 
proposed a basis of master integrals to expand said Green functions, some of them new in 
literature. Due to the potential applications of two- and three-point Green functions with a 
massive gauge boson in physical processes in several frameworks of particle and condensed 
matter physics, we evaluated these new integrals in arbitrary space-time dimensions. The 
basis expansion contains the ordinary expansion of QED3 as the particular case of our 



findings in the parity preserving limit 25|, which is achieved when both the CS term and 



the Haldane mass are taken to zero. Massless theory for both fermions and photons 24 1 
follows as well. Green functions we have obtained will help us make progress in certain 
physical applications in the near future. Particularly, we are interested in the radiative 
generation of fermion masses by the CS term and its impact in several physical observables. 
This project is in progress and the results will be presented elsewhere. 
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Appendix 

In this Appendix we provide details in the calculations for three sets of master integrals: 
for the one- loop photon propagator described in ( fT3|) . the one-loop fermion propagator, 
Eq. f|T8|) . and for the one-loop vertex described in i^^. Let us begin with the -R'-type 
master integrals. 



1. K 



We begin solving K. Using 



dx 



ah Jq [xa + (1 — x)h] 



2 ' 



we can write 



K = I dx I d^w- 



x[{w — pY + + (1 — x){w'^ + m\) 



(25) 



(26) 



The denominator in the above expression can be written as w'^ — 2xw ■ p + xp^ + m^, which 
after the change of variables w' = w — px, allows us to write 



K 



dx I dPw' 



1 



/o J [{w'Y + p'^x{l — x) + m^'^ 
This integral has the following form: 

In our case, we take D = n = 2 and s = —p'^x{l — x) — and then 



(27) 



(28) 



K = n^'^v(^ J\x[-p^x{l 



X) — m 



2. K., 



2 1-1/2 



27r2 / p 



arctan 



p 



2m± 



(29) 



Now, we calculate K^. This integral depends only of p^, thus we can write 



= cp^ . 



(30) 
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We find c multiplying with on both sides 



w ■ p 



(w^ + m^) [{w — pY + 



Using w ■ p = w + p — {w — p) \ we have 



c = -K . 
2 



3. K, 



Finally, we solve K^y. In its most general form, it depends only on p^jPv and 



'-S„ 



p^ 3 

which is symmetrical under the exchange of fi and u. After simplifying, we find that 



^ f J3 1 

a = - / aw 



+ 4m| 



+ 



and 



1 



w 



m 



3 J {w — pY + m\. ' 3 
where the tadpole integrals [see Eq. ( 137|) below] are 

1 f d^w 



d-^w 



w'^ + m\ J {w — pY + 



Next, we provide details of the calculation of the integrals in Eq. fllSp . 



The integral 



4. r{ml) 



d^w- 



1 



can be evaluated directly from Eq. fl28|) with n = 1, D = 3 and s 

r(m|) = -27i'^m±. 



—m\. It yields 
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5. X(^2^m|,p2) 



Let us evaluate 



Using Feynman parametrization, and after the appropriate change of variables, we 

1 




TT dx{x{l - x)p + m^x + 9 {1 - x)) 2 
Jo 



6. X2(m^,p2^ 



In order to evaluate 
we introduce the auxiliary integral 
such that 



w,. 



w'^[{p — wy + m±] 



1 



Using the identity 



1 



r(n)r(p) Jo ' ' [xA + (1 - a;)5]"+P ' 

with n = 1, p = 2, followed by the shift w' = w — px, we can write 

performing the integration, 

vr^ /""^ 1 1 3 

X^(m|,p^) = — / dx{x~^ - x^)\p^{l -x) + m|]"5 , 

^ Jo 

and the remaining integration yields 



X^imi^.p^^ 



rPi^ 



m± 



p2 + m± 



+ I{p,m±) 
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Thus 

X2 m±,p = — = _. 48) 

Next, we will give the results for the new integrals that are described in (12^ . 



7. jW(e,m±) 

We begin solving J^^\9,m±). For this purpose, we consider the following integral in 
Minkowski space: 

Using Schwinger parametrization 

[Ai + ylz]'' 2m J_i^ r(z/) 
to each of the factors in the denominator, the Mellin-Barnes representation for 
JuuU2,u3{6,m±) becomes 



1 



+joo f+ioo f+ioo 



loo J —too -J —too 



xr(-A3)r(-A34 + A3)r(-A5)r(A3 + iyi)r{Xu -h + i^3)r(A5 + 1^2 

t(0) 

'i/l+A3,V2+A5,I'3+A34-A3 



^^ji+A3,V2+A5,I^3+A34-A3(0' 0) ' (^^) 



where we have used the notation A34 = A3 + A4 and dX^ = dX^^ , and where 
"^!^i+A3 U2+X3 i^3+A34-A3(0' 0) been calculated in Ref. 28|. Replacing these results and sim- 
plifying, we obtain the following expression in terms of five integrals to evaluate 

(j\D (u2\D/2-ui23 r r r r /■+*°° 
J^^u.,.Mrn^) = ^3^^r(^.) (27rz)^ J J J J ' d^^dX,dX,dXudX,T{-X,) 



too 



X r(-A.)r(-A3)r(A3 - A34)r(-A5) ( 'Y - ^ ' 



r (A12 + A34 + 1^123 + A5-f)^,, , , D 

^ ^ -r(Ai23 + z^i)r - Ai3 - 1/12 - As + 



Y(D - z/123 - A34 - A5) ^ — - ^ - 2 

xr( - A2 - A34 - ^^13 + ^ ) • (52) 



2 ^ 

We perform one integration at a time. First we integrate over A3, defining 

/a3 = ^ / rfA3r(-A3)r(A3 - A34)r(Ai23 + l^l)T ( - Ai3 - Z/12 - A5 + - ) . (53) 
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Im(u) 



( • ) 



I 

I 

' N I / ■ 

y 
1 



Re(u) 



FIG. 4: Mellin-Barnes contour of integration which separates the poles of r(. . . + u) from those of 
r(. . . — u). We close the contour to the left or to the right, choosing one of the two series of poles. 
Ai and A2 represent the arguments in the sum which are given a Schwinger representation ()50p . 

Using the Barnes lemma 

— / dsT{a + s)T{b + s)T{c - s)T{d - s) = ^ ' ' ' , (54) 

^TTZ J_joo L [a + + c + d) 

with a = — A34 , 6 = z/i + A12 , c = , = — ^^12 — — A5 , we find that 

r(-A34)r(i', + Ai2)r(A2 - i/j - A.5 + f ) , , , 
r(A,-A«-..-A, + f) r(-A,-A„-A,-.,, + -). 

Then, after integrating over A3, Eq. f lS^ becomes 

:7'il,.3(^'"^±) = / L,A4 / / / / C^AidA2rfA34rfA5r(-Ai)r(-A2) 

r(Ai2 + A34 + ^^123 + As - f )r( - A2 - A34 - jyi3 + f ) 

T{D - Z/123 - A34 - A5) 

^ r(-A34)r(^/i + A2)r( - Ai - A34 - A5 - 1/12 + f ) 

r(A2-A34-Z^2-A5 + f) 

xr(A2-z/2-A5 + y) . (55) 

One particular limit of our interest is 6* — )■ 0. We integrate over A5 first and call Ix^ the 
relevant term in the above multidimensional integration. Therefore, a wise selection of the 
position of the poles according to the possible closure of the contour of integration in Fig. H] 
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must be done in order to avoid unnecessary analytic continuation. We close the contour to 
the left, such that | — 6"^ /k'^\ > 1. Then, the poles are due to r(Ai2 + A34 + 1^123 — -f + A5) 
alone, and are located at A5 = — A12 — A34 — z/i23 + y — n . Therefore 

00 i^-lY / £'/2-Ai2-A34-J'i23-nr(A2 + Z^3 + "-)r^Ai2 + A34 + Z/123 - -f + 

" ^""'S^^V r(# + Ai2 + n)r(2A2 + Ai + z/i3 + ^) 



n=0 

Hence, we can write 



(A\D (h.2\D/2-vx2S °° ( — f+ioo f+ioo f+ioo 



ni 
=0 



joo J —too J —too 



. . ,r(A2 + Z/3 + ^)r(Ai2 + Z/l)r(Ai2 + A34 + Z/l23-f + n) 

xi (-A2ji (-A34j — r~r^rr7 — \ T"; — ToT^; 

+ n + Ai2jr(r2 + i^^ + Ai + 2A2) 



Xr(-A2-A34-^/l3 + ^)^py y^y y ^ 

(/)2\ D/2-Ai2-A34-j^i23-n 
- f) • '■'^^ 

Next, we integrate with respect to A34 closing the contour in Fig. H] to the left, where 
I — im? /k'^\ > 1. Naming Ix.^^ the relevant term, we have that 



'A34 



. - T^) ( - 

■400 



rfA34^ - -^j - r(-A34)r(^ - A2 - A34 - ^13 + y 



xT ( A12 + A34 + Z/123 - y + n 



The poles in Ix^^ are due to r(Ai2 + A34 + z/i23 — -2 + ^^^d located at A34 = — A12 
1^123 + -D/2 — n. Thus 

°° ( — ^V / /32\'/ -^2\ -^i2~!^i23+£'/2-n-Z 
Xr|^Ai2 + / + 1^123 

Then we can write 



lOO ^ — JCXD 



xr(-A2)(?, --^ r(Ai2 + .i)r(A2 + z.3 + n) 



xr A12 + z/123 - y + n + / , r,' 1 ^ ^ • (57) 

2 / r(Ai2 + Y + rijr(2A2 + Ai + n + 1/13) 
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Now we integrate with respect to A2, but now closing the contour in Fig. IHto the right, such 
that \(f < 1. Calhng Ix^ to the corresponding contribution, we have 

"+ioo /„2\ -'*2 / ™2\ -\i2-i^i23+D/2-n-l 



'A2 



d^^l^j (-^) r(-A2)r(Ai2 + //i)r(A2 + z/3 + n) 

r(Ai2 + n + / + z/i23 - f) 
r(Ai2 + f + n)r(2A2 + Ai + n + 1/13) ' 

There are poles in this expression are due to r(— A2) and are located at A2 = j- Therefore 

h2 = E ^ (f, j ( - ^ j r(Ai + u,+ mn + J + u,) 

r(Ai+n + /+j + z/i23-f) 



r(Ai + n + j + f)r(Ai + n + 2j + z/i3) ' 

and thus 



rfAir(-Ai)(|^J (-^] T{\, + n + l + U2) 

T^M , ■ , ^ r(Ai+j + z/i23- f + n + 

xTfAi + ? + i/i)^ P7^- — -. (58) 

^ '^r(Ai + n + j + f)r(Ai + n + 2j + z/i3) ^ ' 

Finally, we integrate with respect to Ai closing the contour to the right, such that < 1, 

and observing that the poles of the expression are due to r(— Ai) and they are located at 
Ai = s, we find 

^f-^) f-^) r(z/3 + n + j)r(z/2 + n + / + s)r(i/i+j + s) 

r(f + n + j + s)r(i.i3 + ^ + 2j + s) ■ 

Changing to Pochammer symbols according to 

r(z/3 + n + j) = (i/3;n + j)r(z/3) , 

r(z/2 + n + / + s) = (1/2; ■'T' + ^ + s)r(i/2) , 

r(z^i+j + s) = (i/i; j + s)r(z/i) , 



X 
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r 1 2^123 - — + n + l + J + s 



D 



^^123 - + l + i + s\v{ Z/123 - — 



r(i^i3 + n + 2j + s) = (z/13; n + 2 j + 5)1(1/13) , 



D 



\-2 

we reach at the expression 



q_ 



,2 \ i 



m' 



x(i/i;j + s) 



' (^^123 - + I + s + j){u2;n + I + s){u3;n + j) 
(f;r2 + j + s)(z/i3;n + 2j + s) 

(60) 



or, equivalently 



(,)^(_^2)D/2-.i,3 



rK3-f) ^ 
r(f)r(z.i23) ' 



^^123 - f ,'^2,'^3,i^l 



9^ ^ 

9 ' 9 ' 9 ' 9 

m"^ m'^ 



(61) 



where Tq 



is the generahzed Lauricella function of 4 variables. 



0-1 5 (3,2, 03, 0-4 

&i;&2 

This is the general expression in arbitrary space-time dimensions. In our particular case, we 
take z/i = z/2 = 1/3 = 1 and D = 3, and Wick rotates to Euclidean space. We hence obtain 

r(3-|) 



J(0)(^,m±) = (m2)-3/2 



r(i)r(3)„ 



E ■ 

iJ,s=Q 



n\l\j\s\ \m 



X 



{3 - l;n + I + s + j){l;n + I + s){l;n + + s) 



+ j + s)(2;n + 2j + s) 



(62) 



or 



-111 
^■2 

2' ^ 



^ ^ 

9 ' 9 ' 9 ' 9 

m"^ 



(63) 



This completes the calculation of J^^\9,m±). Next we check some particular limits. 
First we consider 6^ = m?. Then 

^m±j-^mj ^2_.^ ^m^yl (^m^yl l^m^yl r(2n + 2j + 2s + 3) 

r(n + j + l)r(^n + j + s + , (64) 
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that can be written as 



- 1 1 1 






p2 


3 




9 ' 





(65) 



where 



is the generahzed Lauricella function of three variables. 

28] where, in our case, 



ai, 02, 03, 04 

This expression coincides with the Eq. (37) obtained in Ref. 
v\ = V2 = V'i = \ and p, ^ yU, u-^ — z/. 

Another hmit of interest is when 6* = 0. This hmit is achieved by returning to Eq. fl52|) . 
where we take directly ^ = and integrate using the Barnes theorem and integrating in 
accordance with above reasoning. Then, 



j(o)(m±) = (m2)-3/2 ^ 



-iY+'+^ [ k^y ( p'Y ( q^^^ 



n,s,j=Q 



nlsljl 



or 



Tis+j + l)T{n + s + l) 
2jT{n + s + j + l)T{n + s + 2j + 2) 



(66) 



^ 1 1 1-i 

21 -^1 -^1 -^1 2 


P 


p2 




^■2 

2'^ 


9 


9 ' 





where $2 



Xi,X2,X3 



Ol, ^2, fll 

tion. This expression is in agreement with Eq. (32) of the Ref. 
ui = 1^2 = I's = 1 and ui p, 1^2 ^ /U, u. 



(67) 

is the three-parameter generalized Lauricella func- 

after the identification 



Finally, in the limit when 9'^ ^ in the Eq. (l62l) . we obtain 



2\-3/2 



E ■ 

n j,s=0 



'-lY+i+s f k^y f q^y f p^\' 



nlf.sl 



^ nj+s+1 



xrn+j+s+- 



r(n+j + l)r(n + g + l) 
2 ) V{n + ] + s + l)V{n + 2] + s + 2] 



+ 



(68) 
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8. Jl^\e,m±) 



We write J^^\0^m±) in its most general form 



that is symmetrical under the exchange of k and p. Here 
1 



JA{k,p) 



A2 



p\k' + mi- e') - {k ■ p){p' + mi- e') 



Jo(6',m±) 



(69) 



+ ik ■ p) X{e\ mi, k') - p' X{e\ mi,p') + -{p^-k-p) Ko(m±) 



JB{k,p) = JA{p,k) 



(70) 



and Joie,m±) = 2/7r^J^°\e,m±), Ko{m±) = 2 / -k^ K'^^\m±) , i{e^,ml,p^) = 
X(^^,m^,p^)/7r^. The limits 6* — 0, 'm± — )■ and 6 — )■ m±. can be derived straight- 
forwardly and show complete agreement with the findings of Refs. 18|, \24\. 



9. 4V{e,m±) 



In its most general form, we expand 



jS){0,m^) = y| - ^-fMm±) - (^k,k, - Vy) Mk,p) 



2k ■ p^ 

Vi^K + kf,py - 5^,.—— ) Joik^p) 



where 



Jc{k,p) 



A2 



-{k ■ p){p' + - r) - p'(r + mi-e 



JA{k,p) 



-^{p' + ml-e')Mk,p) 



(71) 



+ {k ■ p + p ) 



2 Ko{m±) 
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JD{k,p) 



2A2 

+ 



3 



Ja 
2 



\{p-kf + \{k-p) 



3 

i^o(m±) 



Jb 
2 



JE{k,p) = Jc{p,k) 



and X^(^2, m|,p2) = X^{6'^ , m\, p^) / n"^ , being 



This integral is readily found to be 



X,(^2^^^^p2 



T^'^Pv 



I{p, m± + 



(72) 



(73) 



(74) 



The limits ^ — ?■ 0, m± — >■ and 6 — )■ m± can be derived straightforwardly and show 
complete agreement with the findings of Refs. 18|, \24\- 



10. 4'H^,m±) 



In order to integrate 

we use the identity 
Thus, we can write 



w2[i(;2 + 6*2] [(p — wy + [(/c — m;)2 + rri^] ' 



W^[W^ + ^2] 02 



yj2 yj2 



+ ^2 



1 



L«(^,m±) = -[J«(m±)-J«(^,m±)]. 



(75) 

(76) 

(77) 



11. L£^(^,m±) 



21 



Similarly, we can express 

= hjfK^±) - jfKO.^±)] . (78) 
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